Structural properties of silicon dioxide thin films densified by medium-energy particles 
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Classical molecular-dynamics simulations have been carried out to investigate densification mecha- 
nisms in silicon dioxide thin films deposited on an amorphous silica surface, according to a simplified 
ion-beam assisted deposition (IBAD) scenario. We compare the structures resulting from the depo- 
sition of near-thermal (1 eV) Si02 particles to those obtained with increasing fraction of 30 eV Si02 
particles. Our results show that there is an energy interval — between 12 and 15 eV per condensing 
Si02 unit on average — for which the growth leads to a dense, low-stress amorphous structure, 
in satisfactory agreement with the results of low-energy ion-beam experiments. We also find that 
the crossover between low- and high-density films is associated with a tensile to compressive stress 
transition, and a simultaneous healing of structural defects of the a-Si02 network, namely three- 
and four-fold rings. It is observed, finally, that densification proceeds through significant changes 
at intermediate length scales (4-10 A), leaving essentially unchanged the "building blocks" of the 
network, viz. the Si(Oi/2)4 tetrahedra. This latter result is in qualitative agreement with the mech- 
anism proposed to explain the irreversible densification of amorphous silica recovered from high 
pressures (~ 15-20 GPa). 

PACS numbers: 61.43.Bn, 68.55.Ac, 77.55. -ff, 81.15.Aa 
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I. INTRODUCTION 

In response to the rapidly increasing demand for 
highly-specialized applications in optical filtering devices, 
integrated electronics and photonics technology, consid- 
erable efforts have been devoted in the past few decades 
to the development of methods for growing thin silicon 
dioxide (Si02) coatings having low stress and low defect 
concentrations, and exhibiting bulk-like optical proper- 
ties. To this end, methods commonly referred to as ion- 
assisted deposition (IAD), in which growth is accompa- 
nied by|-£i|-|flux of energetic particles, are of particular 
interestn □ indeed, bombardment with low to medium- 
energy neutral or ionic particles — in the range 10-200 
eV — is found to have beneficial effects on a number 
of physical properties of the deposited material. Be- 
sides fulfilling the above requirements, it also leads to 
reduced moisture absorption, higher density, lower chem- 
ical etch rate, better adhesion to the substrate, and there- 
fore better overall long-term behavior.El^Ll In addition, 
IAD methods allow processing at lower substrate tem- 
peratures, making them especially adapted to growth 
on polymers or graded-index coatings .El While there evi- 
dently exists a threshold energy below which benefits are 
insignificant, too high an energy can have undesirable 
consequences, such as preferential sputtering or impurity 



incorporationJj 

In ion-beam assisted deposition (IBAD) — a special 
case of IAD — an external source provides a well colli- 
mated flux of ions with a narrow kinetic energy distri- 
bution. It has proven useful to discuss the dependance 
of structural changes of IBAD films upon bombardment 
in terms of three parameters: the mean kinetic energy 
of the ions, Ei, the ratio of ion-to-condensing particle 
fiuxes at the surface, 4'i/4>n (where 0„ is measured from 
the growth rate), and the energy deJiijered by each par- 
ticle that condenses, Ep — Ei4>i/(j3naBB It has been ob- 
served that there exist "critical" values of these param- 
eters below which the filrps are porous and above which 
they are densely packed.H In the latter case, films are 
in general found to exhibit better physical properties, 
closer to their bulk, relaxed counterparts. For Si02, Ep, 
determined from optical index, step coverage or optical 
transmittance measurements, is in the range 10-100 eV, 
depending on th^jjiaxticular set of parameters used in the 
growth process. fflu'EI 

Clearly, bombardment profoundly affects the mi- 
crostructure of the deposited material. For amorphous 
solids — which is the case of Si02 films — static dis- 
order hampers the measurement of detailed information 
beyond the second-neighbor shell. As a consequence, 
experimental results must be considered in the light of 
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structural models. In this way, proper conclusions on the 
atomic mechanisms that mediate bombardment-induced 
densiiication can be drawn. Thus, for instance, models 
help understanding how changes in the medium-range 
structure of amomhous silica (a-Si02) influence the Si- 
0-Si bond-angle,tj which in turn is related to such prop- 
erties as the den^tjL the optical index or the amount of 
built-in stressa#LlO 

It is useful to discuss the structure of a-Si02 in terms 
of short-range ordep-|(SRO) and medium-range order 
(MRO) correlations jl3 at length scales < 4 A and in the 
range 4-10 A, respectively. SRO is well defined in terms 
of such quantities as bond lengths and coordination num- 
bers. It is mainly associated to the structural "building 
block" of the networkriiie_Si(Oi/2)4 tetrahedcprL. Both 
numerical simulationsE3nla and experimentll3il3 have 
shown the SRO correlations in permanently-densified 
(20%) a-Si02 recovered from high pressures (~ 15-20 
GPa) to be essentially the same as in the "normal" ma- 
terial. 

In contrast, MRO correlations, which are at the origin 
of the so-called "first sharp diffraction peak" (FS. 
undergo important modifications upon compression:E3 
the height of the FSDP decreases and its. position in- 
creases, with density. Molecular-dynamicsEJ and Monte- 
CarloEj simulations have attributed this effect to mod- 
ifications in the topology of the network; in particular, 
the intensity of the FSDP is reported to be roughly pro- 
portional to the number of six-membered rings. (A n- 
membered ring is defined as a closed loop of n Si-0 
bonds). Further, the x-ray diffraction spectra of porous 
silica thin films deposited at low substrate temperatures, 
using either electron-beam evaporatiopEll or low-pressure 
chemical vapor deposition (LPCVD)J£3 exhibit a strong 
low-angle scattering intensity and a FSDP height smaller 
than that of bulk a-Si02. The latter effect, which dimin- 
ishes upon annealing, has been attributed to differences 
at the MRO level. However, to the best of our knowl- 
edge, no systematic studies of the changes in MRO with 
the kinetic energy of the ions or the ion-current density 
have so far been reported; this constitutes one objective 
of the present work. 

Several computational methods have been employed to, 
assess the role of ion beams during thin film growth.c^ 
Among them, molecular-dynamics (MD) simulations 
have provided useful insights into -the understanding pJ 
ion-induced layer-bjs-layer ecpwthjHJ surface mobility,^ 
pore annihilation^^ strestO and defect formation.E^ 
These studies have focused on the role of collision- 
induced events for obtaining crystalline material. Al- 
though the precise values of the parameters required for 
this purpose are material-dependent, a beam energy of 
the order pi A few tens of eV — in agreement with 
experimeniO'LJ — appears to be adequate, yielding or- 
dered structures on the scale of many interatomic spac- 
ings. These simulations allow collision-induced events to 
be sorted out and relevant deposition parameters to be 



identified, thus providing useful information for better 
control of the growth process. 

The situation is more complicated for the case of amor- 
phous materials growth, as spatial order is short-ranged. 
Only very few MD studies of the structural properties of 
amorphous thin film growth have thus far been reported; 
notable exceptions are Refs. for silicon and Refs. 

|3^ , p3| for diamond-like carbon films. In the latter case, 
bonding is found to be partly sp^ and partly sp^ , that 
is, between graphite and diamond. For the case of Si02 
— for which no such studies so far exist — considering 
the large number of polymorphs of the material,E3 as- 
deposited films are expected to show a large diversity of 
structures as a function of density. 

In view of this, and given the technological importance 
of the material, we have carried out detailed MD simu- 
lations of the growth of silicon dioxide on a-Si02 sub- 
strates. More precisely, growth is modelled within a sim- 
plified IBAD scenario: we examine, at the atomic level, 
how the structure of the films evolves as a function of 
_R, the ratio of medium- (30 eV) to low-energy (1 eV) 
Si02 particle fluxes impinging on the surface. The ra- 
tionale for the model is given in Sec. ||. For the sake 
of simplicity, the energetic ions are taken to be Si02 
"molecules" ; although this is a very crude approxima- 
tion to the real growth process, our model is expected 
to provide much-needed information into the fundamen- 
tal aspects of bombardment-induced densification. We 
find, on average, that densification is related to changes 
in the ring structure, in a way which is similar to the 
case of pressure- induced densification of "normal" (2.2 
g/cvo?) silica glass. The intensity of the FSDP and the 
number of six-membered rings are here found to increase 
with density. We also find that the crossover from low- to 
high-density films is associated with a tensile-compressive 
stress transition, occurring for Ep of about 12-15 eV. It 
is observed, finally, that medium-energy bombardment 
promotes the annihilation of structural defects, namely 
three- and four-niembered rings, which are at the origin 
of the D2 and Dj—defect lines observed in the Raman 
spectrum of silica£3 



II. COMPUTATIONAL DETAILS 

MD consists in integrating the classical equations of 
motion for a set of interacting particles. The interactions 
can be derived from first principles, or expressed in some 
effective, empirical form whose parameters are fitted to 
experimental data. The first-principles approach is very 
accurate, but does not allow large systems to be dealt 
with over "long" times cales. These limitations are not 
so much of a problem in the empirical approach, but accu- 
racy and "transferability" , that is, the ability for a model 
potential to adequately describe a phase to which it was 
not fitted, is an issue. For the problem considered here, 
empirical potentials are unavoidable. The model devel- 
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oped by Nakano et al. has been successfully employed to 
describe the bulk and the surface of discpdesed, phases of 
Si02, viz. amorphous and porous silica J13't3'L3cll which 
are evidently relevant to the present problem. In this 
model, silicon is assumed to be coordinated to four oxy- 
gen atoms, that is, the network is chemically ordered. 
The total potential energy of the system is the sum of 
two- and three-body contributions: 

i<j i,j<k 

where i, j, and k run over all particles. The two-body 
term is written as 

-^4i^.-f, (2, 

where the various terms represent the steric repulsion, 
the Coulomb interaction due to charge transfer, and a 
charge-dipole interaction due to atomic polarizability, re- 
spectively. The three-body term, which describes the co- 
valent bonding, is expressed as 

v}'2(.r^J,rJk) = exp(— ^ + ^ ) 

x[cos(0.,fc)-cos(&^.,)]2 

xe(ro -rjfe)e(ro -ry ), (3) 

where & is the Heaviside step function, and O^jf. is the 
equilibrium value of the (rjijrjk) bond angle. This term 
is calculated for nearest-neighbor Si-O-Si and 0-Si~0 
triplets only. The model was fitted to various physical 
properties of bulk Si02, such as cohesive energy, elastic 
constants and the phonon density of states; the values of 
the parameters can be found in Ref. 

As demonstrated by surface extended x-ray absorption 
fine structure and x-ray photoelectron spectroscopy mea- 
surements, the atomic structure of silica coatings stroncJyj 
depends on the particular deposition process usedicaa 
both silicon and oxygen atoms can be found in differ- 
ent environments, allowing the presence of homoatomic 
bonds. This is particularly true when sub-oxides are 
formed. These effects are not considered in our model; 
thus, the role of structural defects — such as peroxide 
radicals and, more generally, homopolar bonds — on re- 
laxation mechanisms are not taken into account. For 
the case of (nearly) stoichiometric situations, experimen- 
tal data on Si02 produced by plasma enhancedj-chem- 
ical vapor deposition LE|ECVD),eI RF sputteringjlp-fiec- 
ondary ion depositionpj and ion-beam deposition^ in- 
dicate that the basic "building-block" of the network is 
the Si(Oi/2)4 tetrahedron. Thus, chemical ordering is 
expected to dominate in these cases, and the model po- 
tential of Nakano et al. is assumed to correctly describe 
atomic correlations. 



In IBAD, silicon (di-)oxide|-|is evaporated in a low- 
pressure oxygen atmosphere;crcl part of the condensing 
particles therefore have an energy smaller than 0.1 eV. 
Simultaneously, a beam of 0^/0^ ions with energy in 
the range 30-500 eV is directed onto the substrate, in 
order to promote surface relaxation as well as to hamper 
the formation of suboxides due to preferential sputtering 
of oxygen atoms. Although the experimental procedure 
is relatively simple, an atomistic description of it is a 
formidable task, as it requires an environment-dependent 
potential for oxygen ions, neutral gaseous particles, and 
for surface and bulk atoms. For simplicity, because large 
systems have to be dealt with, we describe the growth 
of silica in terms of a single type of condensing parti- 
cle, namely Si02, for both the deposited species and the 
incident energetic ions. Since the clusters are free to dis- 
sociate on the substrate (and actually do so), we expect 
this approximation to be of little consequence on the final 
structure of the deposited film. However, because oxy- 
gen atoms are preferentially sputtered as the deposited 
energy is raised (leading to a 10% deficit in the worst 
cases), Si03 clusters were introduced when necessary so 
as to maintain the stoichiometry at 2:1. Also, the en- 
ergy transferred during a collision varies with the mass 
of the incoming particle, which, in turn, may influence 
the value for which dense films are obtained. Thus, care 
must be taken when comparing these results with exper- 
iment. We also neglect atomic ionisation and electronic 
loss during collisions, a reasonable approximatipii_at low 
energies, especially for the case of insulators. The 
present model is evidently crude, but we expect that the 
generic features of the structure will be relevant to real 
materials, even though the dynamics of the deposition 
process is not described precisely. 

The substrate on which growth takes place consists of 
a 27.59 X 27.59 x 10 amorphous slab containing 747 
atoms, extracted from a bulk (27.59^ A^) sample pre- 
pared using the activationptelaxation technique (ART) 
of Barkema and Mousseau.cZl The substrate is first sub- 
jected to an extensive thermal annealing cycle in order 
to allow the surface to relax. A 3 A-thick layer at the 
bottom is held fixed, in order to mimic a semi-infinite sys- 
tem, but also to prevent the substrate from moving be- 
cause of momentum transfer from incoming atoms. The 
remaining 7 A are coupled to a heat bath at room tem- 
perature, consistent with experiment .B'tl Velocity renor- 
malization was used to maintain the system at constant 
temperature. 

Si02 particles of roughly thermal kinetic energy (E — 
0.02 ~ 0.06 eV) were used in an attempt to first grow 
samples without bombardment. However, this led to un- 
stable structures and to clusters of a dozen atoms desorb- 
ing from the growing film. Much better results were ob- 
tained with pE = 1 eV, an energy similar to that used by 
Hensel et aZ.Ell (2 eV) in their MD study of bombardment- 
induced epitaxial growth of elemental silicon. For the 
high-energy particles, we used 30 eV, aa suggested by 
the experimental results of McNeil et alB As noted pre- 
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viously, important process parameters are expressed in 
term of the ion-to-condensing particle fluxes ratio, (pi/cfin- 
However, this quantity is not known a priori {(pn is de- 
duced from the measurement of the growth rate); thus, 
it is not a adequate parameter for monitoring the sijaula- 
tion procedure. Instead, we used, as other authors,Lil the 
ratio, R, of low-to-medium energy particle fluxes. Low- 
and high-energy particles were chosen at random accord- 
ing to the desired R value, and deposited at normal inci- 
dence onto the substrate. The equations of motion were 
integrated using the Verlet algorithm with an adapta- 
tive time step, so as to adiust to abrupt changes of the 
velocities during collisions M3 The time interval between 
the deposition of two successive particles was chosen in 
such a way as to provide a reasonable compromise be- 
tween computational effort and proper relaxation of the 
surface. In practice, we found it adequate to relax the 
system during 4 and 8 ps after depositing a low- and a 
high-energy particle, respectively, before a new particle 
was introduced. 

As the film grows, "normal" thermal diffusion into the 
substrate is not sufficient to drain away the excess energy 
brought about by incoming particles. Therefore, the en- 
tire sample was submitted to a 1 ps thermalisation phase 
each time before introducing a new particle. Prior to 
this, the temperature was around 350 K and 320 K for 
30 eV and 1 eV Si02 particles, respectively. Since the 
constants of diffusion are small at these temperatures, 
the equilibration procedure is expected not to interfere 
with the growth process in any significant manner. 

Simulations were carried out for several values of the 
parameter R, namely 0, 0.08, 0.22, 0.38, 0.55, 0.73, 1.07, 
1.53, 2 and 2.73. It can be determined a posteriori that 
these correspond to (pi/cpn values of 0, 0.07, 0.18, 0.32, 
0.41, 0.49, 0.62, 0.75, 0.80 and 0.9, respectively. In each 
case, about 1300 atoms (~430 clusters) were deposited. 
Following deposition, the samples were equilibrated dur- 
ing 20 ps at 300 K, after which a 40-ps microcanonical 
run was carried out in order to let the sample reach struc- 
tural equilibrium. The corresponding "true" duration of 
simulated growth is about 3 ns. All calculations were per- 
formed using the program groF, a general-purpose MD 
code for bulk and surfaces developed by one of the au- 
thors (LJL). 



III. RESULTS AND DISCUSSION 

A. Structure of samples grown without bombardment 

In order to assess the role and importance of bombard- 
ment on the properties of the grown material, we first 
discuss the case R = 0, where all incoming particles have 
the same energy, viz. 1 eV. Figure |l] shows a snapshot 
of the system at the end of the run. The material is evi- 
dently highly disordered and porous. In Fig. ^, we plot, 
for this system, the corresponding integrated number of 



atoms from the bottom of the simulation cell; we note 
that it increases approximately linearly with z in the re- 
gion 13-53 A. This portion will henceforth be referred to 
as the "bulk-like" part of the deposit, where structural 
analysis will be performed in order to limit artifacts due 
to the surface and the substrate. The same procedure is 
used for samples grown with bombardment, which leaves 
us with 20 to 40 A of bulk material, about 1000 atoms 
in all cases. 

The overall mass density in the bulk-like region of the 
slab, 1.3±0.1 g/cm^, is about 40% less than the "nor- 
mal" density of a-Si02 (2.2 g/cm'^) and 12% less than 
the density of Si02 grown by electron-beam es^aporation 
on a silicon substrate at room temperature. □ The dis- 
crepancy between experimental and calculated densities 
is possibly due to the short duration of the simulation 
runs, which may not be long enough to allow complete 
relaxation, and thus full densification. (The computa- 
tional effort required to achieve complete relaxation is 
formidable, and it constitues one of the outstanding prob- 
lems of disordered-system simulations). For this model 
of silica, it is worth noting that the density is close to 
the critical value of 1.4 g/cm^, belop, which pore perco- 
lation occurs and leads to fracture.c2l This was not ob- 
served here; however, as will be shown below, the film 
contains high tensile stress (~ 2 GPa). Thus, relaxation 
is expected to proceed by densification and/or fracture 
formation. This is consistent with the observation of re- 
duced crack toughness and adhesion to the substrate for 
the case of coatings grown at low substrate temperature, 
Ts, with insufficient bombardment .□ 

We now discuss the structure and topology of the de- 
posited material, in order to provide a reference for sam- 
ples grown in the presence of bombardment (cf. Sec. 
IIIB). Figure || shows the total pair distribution func- 
tion for the bulk-like part of the sample, as well as for 
bulk silica; also shown are the partial distributions .gsiSi, 
gsio and goo- Here again, the amorphous character of 
the material is clearly evident: order beyond the nearest- 
neighbour shell rapidly fades out. The nearest-neighbor 
Si-0, Si~Si and 0-0 distances-jare hsted in Table |, 
along with those of bulk a-Si02 Jl3 as well as other mate- 
rial produced-by secondary ion deposition (SItti at room 
temperatureN and LPCVD at Tg = 430° C;B the cal- 
culated and experimental values are seen to be in satis- 
factory agreement. The partial pair correlations reveal a 
broadening of the first peak of the Si-Si distribution com- 
pared to the "normal" glass, presumably a consequence 
of the more disordered structure of the deposited film. 
Also, small pre-peaks to the left of the nearest-neighbour 
peak are seen to develop in the Si-Si and 0-0 correla- 
tion functions, at 2.45 and 2.32 A respectively. As will 
be demonstrated below, these are due to the presence of 
two- and three-membered rings, leading to Si-O-Si and 
0-Si-O bond angles of 90 and 130 degrees, respectively. 

In spite of the disordered nature of the network, the 
structural building block remains, to a good approxima- 
tion, the Si(Oi/2)4 tetrahedron, even though it might 
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be slightly distorted. Indeed, the coordination of silicon 
atoms, calculated from the area under the first peak of 
gsio, is about 3.9. Further, as displayed in Fig. ^(a), 
the 0-Si-O bond angle distribution has a major peak 
centred at 109.7°, the usual tetrahedral bond angle. 

The topology of the network can be described more 
precisely in terms of rings of nearest-neighbor bonds, 
where an n-membered ring is defined as the shortest 
closed path, of n Si-0 bonds originating on a given sili- 
con atom.Ej In normal silica glass, the statistics are dom- 
inated by five-, six-, and seven-fold rings£3 Smaller rings 
are unlikely to occur because they are too expensive in 
terms of elastic bond-bending energy. The existence of 
"odd-numbered" rings results in a broad distribution of 
Si-O-Si bond angles. In the normal material, this distri- 
bution is centered around 143° and has a width of about 
25°j£d as can be seen in Fig. §(b). In contrast, the mean 
Si-0~Si bond angle for evaporated Si02 is known from 
infrared spectroscopy measurements to be 7-10° smaller 
than in a-Si02.a This is attributed to the presence of 
three- and four-membered rings, manifested by the pres- 
ence of two sharp peaks at 606 and 495 cm~^ in the Ra- 
man spectrum of chemical vapor deposited (CVD) films, 
produced at low substrate temperature.E3 

Such small rings are a common feature of disordered 
Si02 networks, . and. they are a consequence of rapid 
quenching rates;t§E3 thus, they are also expected to be 
present in Si02 produced using evaporation methods. 
In the present model, the fraction of three- and four- 
menibered rings is about 0.1 and 0.15 per deposited 
silicon atom, respectively, roughly ten and two times 
that of normal a-Si02.Ej The mean value of the Si-0- 
Si bond angle is 139°. These results are consistent with 
the above-mentioned experimental results for evaporated 
films. In addition, two-membered rings (edge-sharing 
tetrahedra) are present here, but absent in the normal 
material. As can be seen in Fig. ^(b), the presence of 
such highly-strained rings gives rise to a sharp peak at 
low angles, ~ 92°, in the Si-O-Si bond angle distribution, 
while the distributions at high angles agree within statis- 
tical accuracy. The 92° peak is directly related to two- 
mcmbered rings, as can be appreciated from Fig. ^(c), 
while three-membered rings give rise to a shoulder at 
about 130° in the bond-angle distribution. 



B. Structure of samples grown unth bombardment 

We now examine the structure of samples produced in 
the presence of bombardment. It is expected that the 
energy provided by energetic particles will have an "an- 
nealing" effect, thereby reducing the proportion of high- 
energy structures — mostly two- and three- fold rings. 

Simulations were carried out for several values of the 
R parameter. As an example, we show in Fig. |5|the final 
configuration for R = 0.73. Comparing this with Fig. 
1 the main role of bombardment — densification — is 



clearly evident: the mass density of this new sample is 
p — 2.16± 0.1 g/cm'^, much higher than the value of 1.3 
g/cm^ for the sample grown without bombardment, and 
close to 2.2 g/cm"^ for normal a-Si02. (The density of a- 
quartz is 2.65 g/cm"^). In Fig. ||(a), we plot p as a function 
of R] it is seen to increase rapidly, reaching a plateau at 
about 2.3 g/cm^ for i? > 0.6 = Re- This "critical" R 
value roughly marks the crossover from a state of tensile 
stress (at small values of R) to a state of compressive 
stress, as can be seen in Fig. mh) where the pressure is 
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plotted vs R. Near i?c, the ratio of ion-to-condensing 
particle fluxes 4'i/4'n ~ 0.4 — 0.5; this yields a critical 
energy, i^c, for obtaining high-density, low-stress films, 
Ec — Ei(l)i/(f>n, on the order of 12-15 eV per condensing 
Si02 unit. 

The calculated value of Ec is in reasonable agreement 
with experiment, considering the relatively crude nature 
of our model and the evident difficulty in measuring this 
quantity in a reliable manner. Indeed, values ha«e been 
reported between 10 and 100 eV.Bu Souche et al.p for in- 
stance, gave 35 eV, based on optical index measurements 
on Si02 prepared by oxygen IBAD at T5 ^ 250°C, for a 
beam energy Ei = 150 eV. McNeil et al.p using a sim- 
ilar deposition process but lower Ei, 30 eV, found the 
optical transmittance to deteriorate with high ion cur- 
rent density, indicating that Ec is smaller than 3|&,eV; 
an estimate in the range 10-20 eV has been given.E3 Fi- 
nally, Al-Bayati et al.jtn using 10 eV Si"'"/0+ ion beams 
to grow Si02 on a silicon substrate at 350°C, found the 
interface to be smooth and free of strain, and the film 
to be of excellent quality, with no evidence of suboxides; 
thus, 10 eV would appear to be an upper bound to Ec- 

The structural changes which take place upon densifi- 
cation can be investigated in terms of short- and medium- 
range order modifications. At short range, Si-0 and 0-0 
nearest-neighbour distances are found to depend very lit- 
tle on R, and they are close to the values observed for th e 
sample grown without bombardment (cf. Section III A). 
We note, however, a slight ('^2%) increase of the Si-Si 
distance, consistent with an increase of the Si-O-Si mean 
bond angle from 139° for i? = to 143° for R = 2.73. 
The Si coordination, between 3.90 and 3.96, shows no 
systematic variation with rising density. 

However, the picture is quite different at intermediate 
length scales. In Fig. ^ we plot, for both R — and 
R = 2.73, the total pair distribution function, g{r), and 
the running coordination number, N{r); the latter, essen- 
tially the integral of the former, gives the mean number 
of atoms within a sphere of radius r centered on an av- 
erage atom, irrespective of the atomic species. It is clear 
that the N(r) curves for the two samples are very similar 
for r < 3 A, confirming the role of the Si(Oi/2)4 tetra- 
hedron as the structural building block (no homoatomic 
bonds are allowed in this model). Above this value, how- 
ever, strong differences develop. Densification, therefore, 
involves structural modifications at or beyond second 
nearest- neighbors . 

Amorphous silica is known to possess significant or- 
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der at medium range, manifested in the so-called "first 
sharp diffraction peak" (FSDP) at 1.55 A~^, and oh-, 
served in neutron and x-ray diffraction experiments. E2I 
Although the precise origin of this peak,|iH|terms of sta- 
ble entities, is still the matter of debate,E3"E3 MD simu- 
lations by Nakano et al. have shown that its position and 
shape arise from_cprrelations between atoms at distances 
as large as 12 A.E3 The FSDP can therefore be viewed as 
a signature of the presence of MRO in a-Si02. In order 
to assess this picture within the present approach, that 
is, to quantify possible modifications in MRO upon den- 
sification, we calculated the static structure factor, S{q), 
which is also accessible experimentally. This is shown 
in Fig. ^ for the samples with R = and R — 2. 23, as 
well as for a-Si02 prepared using the ART method;cZl the 
curves were smoothed using Gaussian functions of width 
0.1 so as to mimic thermal as well as experimental 
broadening. There is evidently very little difference from 
one sample to another for q- values larger than 4 A^^, 
consistent with the existence of a structural unit com- 
mon to all three networks. In contrast, the low-g portion 
of the spectra exhibit significant differences. In particu- 
lar, for q ^ 1 A~^, the intensity varies strongly with _R, 
in a manner which is inversely related to the density. 

While the finite size of the system, plus periodic 
boundary conditions, artificially contribute to the inten- 
sity at such low-q values, they cannot be solely respon- 
sible for the observed variations. Since all three samples 
are treated in the same manner and are similar in size, 
these variations must be related to the presence of pores 
in the low-density samples, an effect which can be under- 
stood in term of a Babinet argument — the scattering by 
a pore is similar to that from a cluster having comparable 
dimensions. 

Changes in the MRO upon densification are also ap- 
parent in the FSDP, which the present model places at 
about 1.60 A~^, in satisfactory agreement with the ex- 
perimental value, 1.55 A~^. For R = 0, that is, highly 
porous material, the FSDP shows up as a mere shoul- 
der, whose intensity is 23% smaller than that of a-Si02. 
For R = 2.73, a well-defined peak is observed, which is 
only 7% smaller than that of the normal material. The 
complete picture is provided in Fig. H, where we plot the 
intensity of the FSDP as a function of density. Although 
the statistics are not perfect, there is a clear tendency 
for the height of the FSDP to increase with increasing 
density. It is interesting to note that if the densification 
were the consequence of a homogenous elastic compres- 
sion of the network, that is, without structural modifi- 
cations, the mean atomic distances would decrease with 
increasing density, implying an overall shift of S{q) to- 
wards higher q values. Here, although the FSDP position 
is modified, no obvious trend is found. 

The differences in heights of the FSDP among the var- 
ious samples are consistent with the results of Himmel 
et aL,tij who found a 28% smaller height in the x-ray 
static structure factor of electron-beam evaporated sil- 
ica, compared with bulk fused silica. Our calculations 



are also in qualitative agreement with MD simulations 
of pressure-induced densification of silica nanophases, 
where the height of the FSDP was found to be 15 % 
smaller coinpared to a-Si02 in the density range 1.67 to 
2.03 g/cm3.El 

The ring structure of the various models can tell us 
more about changes in MRO upon densificationij-because 
the size of rings is typically in the range 4-10 AM modi- 
fications in the MRO are related to those in the arrange- 
ments of neighboring Si(Oi/2)4 tetrahedra. The occur- 
rence of the various types of rings as a function of R, 
relative to i? = 0, is shown in Fig. |o|. In spite of statisti- 
cal uncertainties, the following clear trends can be identi- 
fied, (i) The number of two-, three- and four-membered 
rings decreases with rising R; indeed, as noted in Sec. 
Ill A, these correspond to high-energy structures and 



therefore are expected to anneal out during low-energy 
bombardment, (ii) The number of five-, six-, and seven- 
menibered rings increases with rising R; these rings being 
close to the most probable size in a-Si02, they are thus 
strongly favored. (Hi) The number of nine-membered 
rings decreases with density, indicating a reduction in 
the number of pores. In contrast, the relative occurrence 
of eight-membered rings exhibits no clear trend, most 
likely because they are relatively stable entities in dense, 
cristalline forms of silica, j-such as a-quartz (2.65 g/cm"^) 
and coesite (2.92 g/cm^)^ 

MD simulations on elemental silicon growthEjlEl have 
demonstrated, in agreement with experiments, that low- 
energy ion bombardment has an effect similar to anneal- 
ing. For the case of Si02 prepared by IB AD, our results 
indicate that the elimination of two-, three- and four-fold 
rings may be achieved with the use of a second particle 
beam with kinetic energy in the range 12-15 eV. The con- 
centration of three- and four-fold rings was repcpted to 
decrease for CVD Si02 submitted to annealing;E3 how- 
ever, to the best of our knowledge, no such investiga- 
tion has been reported regarding their dependence upon 
IBAD ion-current density or ion kinetic energy. Such 
data would be of utmost interest. 

The modifications in the ring statistics were pointed 
out in MD and Monte-Carlo sirmilalions of pressure- 
induced densification of bulk silica.lijli3 It was found that 
most topological modifications are mediated by defect 
annihilation and bond switching. The former mecha- 
nism proceeds through the formation of bonds between 
non-bridging oxygen and three-fold coordinated silicon 
atoms, while the latter involves the diffusion of coordina- 
tion defects (hopping of a dangling bond from one atom 
to a neighbor, accompanied by the formation of a new 
bond)lla This has also been identified as an important 
relaxation mechanism in a-Si02, and it appjears to be 
characteristic of chemically-ordered networks O The lat- 
ter point being an assumption of our model, it is expected 
that these mechanisms also play an important role dur- 
ing bombardment-induced densification. More precisely, 
since coordination numbers do not show systematic varia- 
tions with increasing density, bond-switching is expected 
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to dominate over the annihilation of coordination defects. 

In pressure-densified material, the evolution of the ring 
populations is quite different from the one we observe: 
the number of four-, eight-, and nine-fold rings increase, 
while those of five- and six-fold rings, as well as the inten- 
sity of the FSDP, decrease with increasing density (in the 
range 2.2-3.2 g/cm'^). Taken together with the present 
results, this may indicate that in its most stable form, 
a-Si02 has maximum MRO at ambient temperature and 
pressure. 

IV. CONCLUDING REMARKS 

In this article, we have presented the results of ex- 
tensive and detailed molecular-dynamics simulations of 
the growth of Si02 on a a-Si02 substrate, with a view 
of understanding the process of densification induced by 
low-energy particle bombardment. The physical picture 
which emerges from our observations is the following: the 
network undergoes significant density variations through 
structural changes at an intermediate length scale; more 
precisely, ring statistics are affected, while the structural 
unit of the network — the Si(Oi/2)4 tetrahedron — re- 
mains essentially unchanged. Upon bombardment, small 
rings, which are associated with defects in the network, 
rearrange to produce larger ones, thus reducing the en- 
ergy cost associated with bond bending and increasing 
the mean Si-0~Si bond angle. Large rings, particularly 
nine-membered ones, "disintegrate" into five-, six-, and 
seven-membered rings, thus reducing the occupied vol- 
ume and causing the density to increase. Simultaneously, 
this influences the intensity of the first sharp diffraction 
peak (FSDP), which is found to increase with rising den- 
sity. By calculating the pressure inside the sample, we 
have found an "optimal" incident energy range for low- 
stress, high-density films of 12 to 15 eV per condensing 
Si02 particle. 

The very brief (~4 ns) real-time equivalent duration of 
our simulations may not be sufficient to allow the samples 
to escape from local minima in the potential energy sur- 
face, leading to structures which are incompletely relaxed 
compared to the real material. Because relaxation pro- 
cesses, which are expected to increase MRO correlations, 
may occur within the time scale of an experiment, control 
over the film microstructure may be a difficult objective 
to attain when performing ion-assisted deposition exper- 
iments. In addition, most of MRO modifications arise 
from low-energy bombardment. For this_reason, and as 
demonstrated by experimental results,Bc3 we emphasize 
that any deposition process where the average ion kinetic 
energy is greater than about ten eV should be avoided. 
Besides ion-beam experiments, it was recently demon- 
strated that magnetron sputtering processes could pro- 
duce low-energy high ion flux (up to 30 ions per neutral 
condensing particle), allowing amorphous silicon fllms 
with a tunable amount of MRO to be grown.Ell We ex- 



pect that such findings will stimulate futher experiments 
on Si02 films, which will allow our results to be assessed. 
For example, a systematic study of the peaks associated 
with three- and four-membered rings in the Raman spec- 
tra of silica coatings, and the evolution of the FSDP, with 
the kinetic energy and flux of ions could help enhance our 
understanding of bombardment-induced densification. 

Acknowledgments ~ A. L. is grateful to Gilles Dennler 
for useful discussions. Normand Mousseau is acknowl- 
edged for providing the ART model. This work is sup- 
ported by grants from the Natural Sciences and Engi- 
neering Research Council of Canada (NSERC) and by 
the "Fonds pour la formation de chercheurs et I'aide a 
la recherche" (FCAR) of the Province of Quebec. We 
are indebted to the "Reseau quebecois de calcul de haute 
performance" (RQCHP) for generous allocations of com- 
puter resources. 



* To whom correspondance should be addressed; electronic 
mail: Laurent. Lewis@UMontreal.C A 
^ L. Martinu and D. Poitras, J. Vac. Sci. Tech. A 18, 2627 
(2000). 

^ H. K. Pulker, Surf. Coat. Technol. 112, 250 (1999). 

^ Handbook of plasma processing technology, edited by S. M. 
Rossnagel, J. J. Cuomo and W. D. Westwood, (Noyes Pub- 
lications, Park Ridge, 1990). 

" F. A. Smidt, Int. Mat. Rev. 35, 61 (1990). 

^ A. Brunet-Bruneau, D. Souche, S. Fisson, V. Nguyen Van, 
G. Vuye, F. Abeles, and J. Rivory J. Vac. Sci. Technol. A 
16, 2281 (1998); A. Brunet-Bruneau, J. Rivory, B. Rafin, J. 
Y. Robic, and P. Chaton, J. Appl. Phys. 82, 1330 (1997). 

^' M. V. Bazylenko, M. Gross, and D. Moss, J. Appl. Phys. 
81, 7497 (1997). 

^ D. Souche, A. Brunet-Bruneau, S. Fisson, V. Nguyen Van, 
G. Vuye, F. Abeles, and J. Rivory, Thin Solid Films 313- 
314, 676 (1998). 

* L. Martinu, J. E. Klemberg-Sapieha, O. M. Kiittel, A. 
Raveh, and M. R. Wertheimer, J. Vac. Sci. Technol. A 
12, 1360 (1994). 

^ J. R. McNeil, A. C. Barron, S. R. Wilson, and W. C. Her- 
mann, Jr., Appl. Opt. 23, 552 (1984). 
J. Rino, I. Ebbsjo, R. Kalia, A. Nakano, and P. Vashishta, 
Phys. Rev. B 47, 3053 (1993). 

G. Lucovsky, J. Manitini, J. K. Srivastava, and E. A. Irene, 
J. Vac. Sci. Technol. B 5, 530 (1987). 

P. Vashishta, R. Kalia, J. Rino, and I. Ebbsjo, Phys. Rev. 
B 41, 12197 (1990). 

J. Rino, G. Guttierrez, I. Ebbsjo, R. Kalia, and P. 
Vashishta, in Materials Theory, Simulations, and Paralell 
Algorithms, edited by E. Kaxiras, J. Joannopoulos, P. 
Vashishta, and R. K. Kalia (Mat. Res. Soc. Symp. Proc. 
408, Philadelphia, PA, 1996), 333. 

J. S. Tse, D. D. Klug, and Y. Lepage, Phys. Rev. B 46, 



7 



5933 (1992). 

L. Stixrude and M. S. T. Bukowinski, Phys. Rev. B 44, 
2523 (1991). 

^® R. A. B. Devine, R. Dupree, I. Farnan, and J. J. Capponi, 

Phys. Rev. B 35, 2560 (1987). 
^'^ R. A. B. Devine and J. Arndt, Phys. Rev. B 35, 9376 

(1987). 

S. Susman, K. J. Volin, D. L. Price, M. Grimsditch, J. P. 
Rino, R. K. Kalia, P. Vashishta, G. Gwanmesia, Y. Wang, 
and R. C. Libermann, Phys. Rev. B 43, 1194 (1991). 
" C. Z. Tan and J. Arndt, J. Non-Cryst. Solids 249, 47-50 
(1998). 

S. Sugai and A. Onodera, Phys. Rev. Lett. 77, 4210 (1996). 

B. Himmel, T. Gerber, and H. G. Neumann, Phys. Stat. 
Sol. A 88, K127 (1985). 

J. Konnert, P. d'Antonio, M. Huffman, and A. Navrotsky, 
J. Am. Ceram. Soc. 70, 192 (1987). 

^■^ Thin Films, Modeling of Film Deposition for Microelec- 
tronic Applications, vol. 22, edited by S. Rossnagel and A. 
Ulman, ( Academic Press, New- York, 1996). 

^"^ J. Jacobsen, B. H. Cooper, and J. P. Sethna, Phys. Rev. 
B 58, 15847 (1998). 

T. Ohira, O. Ukai, T. Adachi, Y. Takeushi, and M. Murata, 
Phys. Rev. B 52, 8283 (1995). 

L. Dong, R. W. Smith, and D. J. Srolovitz, J. Appl. Phys. 
80, 5682 (1996). 

C. -C. Fang, V. Prasad, R. V. Joshi, F. Jones, and J. J. 
Hsieh, in Ref. || 117. 

J. Kitabatake and J. E. Greene, J. Appl. Phys. 73, 3183 
(1993). 

G. K. Hubler, Proceedings of the conference on Beam-Solids 
Interactions for Materials Synthesis and Characterization, 
edited by D. C. Jacobson, D. E. Luzzi, T. F. Heinz, and 
M. Iwaki (Mat. Res. Soc. 354, Pittsburg, 1995), 45. 

B. Strickland and C. Roland, Phys. Rev. B 51, 5061 
(1995). 

H. Hensel and H. M. Urbassek, Phys. Rev. B 58, 2050 

(1998) . 

H. P. Kaukonen and R. M. Nieminen, Phys. Rev. Lett. 68, 
620 (1992). 

H. U. Jager and K. Albe, J. Appl. Phys. 88, 1129 (2000). 
F. Liebau, in The physics and technology of amorphous 
Si02, edited by R. A. B. Devine (Plenum, New- York, 1988), 
15. 

R. A. B. Devine and M. Marchand, Appl. Phys. Lett. 63, 
619 (1993). 

T. Campbell, R. Kalia, A. Nakano, F. Shimojo, K. Tsuruta, 
P. Vashishta, and S. Ogata, Phys. Rev. Lett. 82, 4018 

(1999) . 

P. Vashishta, R. Kalia, A. Nakano, W. Li, and I. Ebbsjo, 
in Amorphous Insulators and Semiconductors, edited by M. 
F. Thorpe and M. I. Mitkova, (Kluwer, Dordrecht, 1997), 
151. 

A. Nakano, L. Bi, R. Kaha, and P. Vashishta, Phys. Rev. 
B 49, 9441 (1994). 

A. Nakano, R. Kaha, and P. Vashishta, J. Non-Cryst. 
Solids 171, 157 (1994). 

A. Nakano, R. Kalia, and P. Vashishta, Phys. Rev. Lett. 
73, 2336 (1994). 
*^ M. 1. Trioni, A. Bongiorno, and L. Colombo, J. Non-Cryst. 



Solids 220, 164 (1997). 

A. J. Forty, M. Kerkar, J. Robison, S. M. El-Mashri, and 

E. N. Farrabaugh, J. Phys. C8(47), 857 (1986). 
"3 S. K. Ray, C. K. Maiti, S. K. Lahiri, and N. B. Chakrabarti, 

J. Vac. Sci. Technol. B 10, 1139 (1992). 
*■* C. F. George and P. D'Antonio, J. Non-Cryst. Solids 34, 

323 (1979). 

A. H. Al-Bayati, S. S. Todorov, K. J. Boyd, D. Marton, J. 
W. Rabalais and J. Kulik, J. Vac. Sci. Technol. B 13, 1639 

(1995) . 

M. Hedstrom and H. P. Cheng, Phys. Rev. B 59, 10701 
(1999). 

N. Mousseau, J. Chem. Phys. 112,960 (2000). 
R. Smith, M. Jakas, D. Ashworth, B. Oven, M. Bowyer, 
1. Chakarov, and R. Webb, in Atomic and ion collisions 
in solids and at surfaces, edited by R. Smith (Cambridge 
University Press, Cambridge, 1997), 247. 
"3 P. A. V. Johnson, A. C. Wright, and R. N. Sinclair, J. 
Non-Cryst. Solids 58, 109 (1983). 

S. L. Chan, L. F. Gladden, and S. R. Elliott, in Ref. ||, 
83. 

A. Pasquarello and R. Car, Phys. Rev. Lett. 80, 5145 
(1998). 

J. M. E. Harper, J. J. Cuomo, R. J. Gambino, and H. R. 
Kaufman, Nucl. Instrum. Meth. Phys. Res. B7/8, 886 
(1985). 

S. C. Moss and D. L. Price, in Physics of disordered mate- 
rials, edited by D. Adler, H. Fritzsche, and S. R. Ovshinsky 
(Plenum, New- York, 1985), 77. 
5* S. R. Elliott, J. Phys.: Cond. Mat. 4, 7661 (1992). 

P. H. Gaskell and D. J. Wallis, Phys. Rev. Lett. 76, 66 

(1996) . 

5^ P. McMillan, in Ref. g, 63. 

" J.E. Gerbi, P.M. Voyles, M.M.J. Treacy, J.M. Gibson, and 
J.R. Abelson, unpublished 



TABLE 1. Nearest-neighbor distances and Si-O-Si mean 
bond angle, calculated for a sample grown without bombard- 
ment, along with those of Si02 produced by secondary ion 
deposition (SID), electron beam evaporation, low-pressure 
chemical vapor depostion (LPCVD), and bulk amorphous sil- 
ica. Bond lengths are given in A; in parentheses are the stan- 
dard deviations. 



Si-Si 



Si-0 



O-O Si-O-Si 



,0 



This work 
SID a-Si02' 
e-beam a-Si02l 
LPCVPjpf iO: 
a-Si02No 



3.0 (0.17) 1.61 (0.05) 2.64 (0.1) 

3.06 (0.02) 1.61 (0.01) 2.64 (0.02) 

3.08 (0.13) 1.61 (0.06) 2.60 (0.1) 

3.08 (0.1) 1.61 (0.05) 2.63 (0.08) 



139° 
144.6° 
136° 

143° 
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FIG. 1. Snapshot of the sample grown with 1 eV clusters. 
The total number of atoms is 2040. Silicon atoms are in light 
grey while oxygen atoms are dark grey. Darker atoms at the 
bottom of the cell indicate the fixed layers in the substrate 
region. 




FIG. 2. Total (integrated) number of atoms as a function 
of the distance from the bottom of the cell for the sample 
grown without bombardment {R — 0). The dashed curve is 
the best linear fit to the data in the bulk-like region of the 
sample (as indicated). 




FIG. 3. Total and partial pair distribution functions cal- 
culated for a sample grown without bombardment {R = 0). 
In the top panel, the dashed line shows the total pair dis- 
tributisn function for bulk amorphous silica prepared using 
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FIG. 5. Snapshot of the samples grown with i?=0.73. The 
total number of atoms is 2040. The effect of bombardment is 
clearly visible when comparing to R=0. 
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FIG. 4. Bond-angle distributions for a sample grown with- 
out bombardment: (a) O-Si-0; (b) Si-O-Si; also shown is 
the corresponding distribution for bulk-amorphous silica pro- 
duced with ART (dashed line); (c) Si-O-Si for the special 
cases of two- and three-membered rings. 
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FIG. 7. Total pair distribution function, g{r), and running 
coordination number, N{r), for the samples grown with R = 
and R = 2.73, as indicated. The inset is a magnification of 
N{r) in the region 1.4-3.5 A. 
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FIG. 6. (a) Density and (b) pressure versus R in the bulk 
part of the various samples; in (a), the dotted line indicates 
the density of bulk silica (2.2 g/cm^); the error bars corre- 
spond to the calculated standard deviations. 
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FIG. 8. Total static structure factor for the samples with 
R = 0, R = 2.73, and for amorphous silica. Data are trun- 
cated at 0.7 due to the finite-size of the simulation cell. 
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FIG. 9. Height of the first sharp diffraction peak as a 
function of the density of the samples. The dashed lipe is the 
corresponding value for a-Si02 prepared using ART.c3 
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FIG. 10. Relative variation of the number of n-membered 
rings per silicon atom as a function of R. 
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